Suppression of superconductivity in high-T c cuprates due to nonmagnetic impurities: 
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We studied the effects of nonmagnetic impurities on high-temperature superconductors by solving 
the Bogoliubov-de Gennes equations on a two-dimensional lattice via exact diagonalization tech- 
nique in a fully self-consistent way. We found that s-wave order parameter is almost unaffected 
by impurities at low concentrations while d x 2_ y 2-wa,ve order parameter exhibits a strong linear de- 
crease with impurity concentration. We evaluated the critical impurity concentration n\ at which 
superconductivity ceases to be 0.1 which is in good agreement with experimental values. We also 
investigated how the orthorhombic nature of the crystal structure affects the suppression of su- 
perconductivity and found that anisotropy induces an additional s-wave component. Our results 
support d ;c 2_ H 2-wave symmetry for tetragonal and s + d x 2_ y 2-~w&ve symmetry for orthorhombic 
structure. 
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Symmetry of the order parameter (OP) in high-T c 
cuprates is important both for understanding the mech- 
anism of superconductivity and also for technological 
applicationsu. For example, d-wave symmetric OP effec- 
tively refutes phonon mechanism and for a device made 
of a d-w&ve superconductor having no gap in energy 
spectrum, no refinement would get rid of the dissipa- 
tion at low frequencies, even at low temperatures. 3d 
metal (Zn, Ni, Al, Ga, Fe, ...) atom substitution for 
Cu atoms in high-T c cuprates may identify symmetry of 
the OP0. It is a well known fact that for conventional 
superconductors having isotropic order parameter, non- 
magnetic impurities with small concentrations have no 
effect on critical temperatureoQ while magnetic impuri- 
ties act as strong pair breakers, and as a result of this 
superconductivity is suppressed very rapidlyLrEI. On the 
other hand, nonmagnetic imparities are very effective 
in anisotropic superconductorsEflij. For a pup. [ super- 
conductor, anisotropy leads to increase in TjBEM and 
the critical temperature suppression rate with increasing 
impurity concentration is proportional to the strength 
of anisotropyE3113. Ualike the conventional supercon- 
ductors, in hole-dopedEj high-T c cuprates both magnetic 
(Ni) and nonmagnetic (Zn) impurities suppress T c very 
effectively. 

Dependence of the superconducting properties (criti- 
cal temperature, order parameter, density of states, ...) 
on impurity or point defect concentration is a subject of 
ongoing research. So far, most of the experiments have 
been performed to investigate effects-of .Zjl and Ni substi- 
tution in YBa2Cu307_5 compoundaljo. Alternatively, 
disorder can also he introduced by creating defects with 
ion irradiatiorOL.il, but in this case affected region is 
often uncontrollable. 

In spite of the complexity of the high-T c cuprates 
(boundaries, defects, ...), insufficient control of the actual 
impurity or point defect concentration, solubility, and ho- 
mogeneity of the distribution of the dopants which may 



TABLE I. The critical temperature T c o and the initial drop 
X = [(Tco - T c (x))/T c0 ]/x in various Zn doped YBCO com- 
pounds. In constructing the table, we used T c (x) values at 
x < 0.04 for which \ i s almost x independent. 
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lead to contradictory data, we can summarize some of 
the experimental results as follows: 

• For YBa2(Cui_ 3; Zn a ;)307_ ( 5 compounds, at small 
concentrations, x < 0.04, Zn ions occupy prefer- 
ably Cu-sites in the Cu02 planes, however for 
x > 0.04 the^i^stitaicnt starts to occupy Cu-sites 
in the chaiiOEilrE^niZl. Since this compound has 
two planes and one chain in a unit cell, for x < 0.04 
the actual (effective) impurity concentration m, i. 
e. the number of impurities per unit cell per Cu02 
plane, becomes 3x/2. 

• The critical temperature decreases linearly with jp-. 
creasing impurity concentration in substitutior£3'E£l 
(for x > 0.04 the drop rate decreases due to partial 
occupancy of Zn at chap- sites) and point defect 
concentration irradiationa experiments. 

• Due to orthorhombicity of YBCO material, Cu02 
planes exhibit an anisotropic behavippand admix- 
ture of d- and s-wave is possibleE3c3~EJ. 

• Zn substitution does not alter the carrier concen- 
tration in Cu02 planes23. 
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On the theoretical side, the existing pair-breaking 
models overestimate the suppression of critical tempera- 
ture and predict an incxeasing slope with increasing im- 
purity concentratior£Ja which contradicts the observed 
linear dependence of T c on n^. The effects of neither 
magnetic (Ni) nor nonmagnetic (Zn) impurities on the 
superconducting properties of the cuprates have been 
explained clearly. Main reason for the discrepancy be- 
tween theory and experiments is that the conventional 
Abrikosov-Gor'kov (AG) -type pair-breaking models ig- 
nore position dependence of the order parameter, near 
impurity sites. Recently,_Eranz and his coworkers^, and 
Zhitomirsky and WalkerE3 argued that spatial variation 
of the order parameter must be taken into account for 
short coherence length superconductors. 

In the present paper, we investigate effects of nonmag- 
netic impurities on high-T c cuprates for both tetrago- 
nal and orthorhombic phasesJpK solving the Bogoliubov- 
de Gennes (BdG) equationstZIH in a fully self-consistent 
way. In particular, we address possibility of extracting 
OP symmetry by examining effects of nonmagnetic impu- 
rities. Our results support d x 2_y2-wave pairing symme- 
try for tetragonal and s + d x 2_ y 2-w&ve for orthorhom- 
bic structure. Possibility of admixture of s- and d- 
wave symmetries, have already been propjosed-jn various 
experimentallj'EJ and theoretical worksE3~EllO. 

The BdG equations on two-dimensional lattice have 
the following formE3'E3 
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where u n (i) and v n (i) are quasiparticle amplitudes at site 
i with eigenvalue E n , and Ay is the pairing potential. 
The normal-state part of the Hamiltonian can be written 
as 



= (t tj + Uijriij /2)(1 - 5ij) 
+(^ imp - /i + U H na/2)Si. 



(2) 



where is the hopping amplitude, fj, is the chemical po- 
tential, UijUij/2 and Uunu/2 are the Hartree-Fock po- 
tentials with on-site interaction Ua and off-site interac- 
tion Ujj, respectively. Finally V^ lmp is the impurity po- 



tential. The pairing potentials are defined by 

ij — UijFij . 



(3) 



The charge density in the Hartree-Fock potentials 
and the anomalous density Fij in the pairing potential 
are determined from 



F, 



(*T(i)*iO')> 



(4) 



(5) 



where a is spin index, and and ^f a {i) are related 

to the quasiparticle creation (7^) and annihilation ("fna) 
operators 
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FIG. 1. Normalized s- and d-wave order parameters, 
(Ad) /Ado and (A s )/A s o, versus impurity concentration m 
for tetragonal structure. Ado and A s o are the magnitudes of 
the order parameters in the absence of the impurities, and 
(• ■ ■) is taken over 20 different impurity ditributions. Solid 
lines represent the best linear fit to the data. 
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where 7 and 7^ satisfy the Fermi commutation relations. 
The self-consistency conditions can be written in terms 
of u n , v n , and E n 



n 

+v n (i)v*(j)[l-f(E n )} , 

F tj =J2u n (iK(m -f(En)] 
n 

-vt{i)u n {j)f{E n ) , 



(7) 



(8) 



where f{E n ) — l/[exp (En/ksT) — 1] is the Fermi distri- 
bution function. 

We solve the BdG equations on a 20x20 square lattice 
(hence, we diagonalize a 1600x1600 matrix) with periodic 
boundary conditions by exact diagonalization technique 
using IMSL subroutines. After choosing a suitable ini- 
tial guess for OP, we solve Eq. (|l|). Next, we calculate 
the new charge density n,j and anomalous density F^ 
via Eqs. (Q) and (||) and iterate this procedure until a 
reasonable convergence is achieved. The BdG equations 
are solved self-consistently. Self-consistency conditions 
[Eqs. (0) and (||)] lead to 10 separate equations. The first 
five [obtained from Eq. (0)] renormalize on-site energies 
and hopping matrix elements while the last five [obtained 
from Eq. (g)] affect the on-site and nearest-neighbor in- 
teraction terms. Although, the first five of these self- 
consistency conditions can be neglected for conventional 
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FIG. 2. Normalized order parameter versus impurity con- 
centration for orthorhombic structure. Ad^o is x component 
of the order parameter in the absence of impurities. Solid line 
represents the best linear fit to the data. 
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FIG. 3. Normalized order parameter versus impurity con- 
centration for orthorhombic structure. Ad y o is y component 
of the order parameter in the absence of impurities. Solid line 
represents the best linear fit to the data. 



superconductors where U/t <C 1, for strong interaction 
case we should keep them, since they play important role 
especially in the presence of impurities. The impurity po- 
tential V- mp is treated in the unitary limit, i.e. V" np » t, 
and taken nonzero for randomly chosen lattice sites. 

We first solve the BdG equations for tetragonal case 
t x = t y = t where t x and t y are nearest-neighbor hop- 
ping amplitudes along x and y directions, respectively. 
For s-wave OP symmetry we assume that on-site (attrac- 
tive) interaction is Uu — — 1.7t and there is no nearest- 
neighbor interaction. In the case of d-wave OP symmetry 
on-site (repulsive) interaction is Uu — lAt and nearest- 
neighbor (attractive) interaction is = — 1.4i. With 
this choice of parameters we fix chemical potential \x so 
that band filling factor (n) ~ 0.8 and zero temperature 
coherence length £ — 4a. These values are in good agree- 
ment with the commonly accepted experimental values. 

Figure 1 shows that, at low impurity concentrations s- 
wave OP symmetry is almost unaffected by the impurities 
or point defects. Although we use several on-site inter- 
action values by keeping the band filling factor and zero 
temperature coherence length constant, we don't get any 
qualitative change. This result, is consistent with Ander- 
son theoremEl and AG theoryB. However, experimental 
data for high-T c cuprates exhibit a much stronger sup- 
pression of superconductivity with increasing disorder. 

On the other hand, our d-wave calculations give re- 
sults similar to the behavior observed in experiments. 
For d-wave symmetry we find that on-site pairing po- 
tential is negligibly small. In Fig. 1, is amplitude 
of the nearest-neighbor pairing potential. We obtain a 
linear decrease in the mean OP, which is assumed to be 
proportional to the critical temperature T C E3, and the 



slope of the straight line is in good agreement with the 
experimental data summarized in Table I. The critical 
impurity or point defect concentration n% at which su- 
perconductivity ceases is also near to experimental value 
~ 0.1. In comparing our results with experimental data 
we should keep the following point in our mind. For 
x < 0.04, substiijutipaal impurities go to preferentially 
Cu02 planesoHcJcZl and hence the actual concentra- 
tion is 3x/2. However for higher concentrations some of 
the Zn atoms occupy the chain sites, and in this case 
we can not relate the in plane concentration to the ac- 
tual one. Therefore, we used the initial points, i. e. 
x < 0.04, to evaluate the initial drop in Table I. To obtain 
the experimental value for critical impurity concentration 
n^, we extrapolated the linear parts of the experimental 
curves to intersect the impurity concentration axes. 

Similar equations have already been solved by Xiang 
and WheatleyEj, however our additional self-consistency 
conditions and choice of parameters lead to a correct pre- 
diction for the critical impurity concentration. It is im- 
portant to note that we can not have a self-consistent so- 
lution of the OP for extended s-wave by using any physi- 
cal values for the above model parameters. This fact was 
pointed out by Wang and MacDonaldE.1, and they found 
that extended s-wave component is smaller than d-wave 
component by about two orders of magnitude. 

When we introduce an orthorhombic distortion by-tak- 
ing t y = 1.5t x , as suggested by experimental dataOca, 
we observe that an s-wave component (of approximately 
ten percent of d-wave components) is induced. Figures 
2 and 3 show the variation of Ad x and A^y components, 
respectively. In the absence of disorder, Ad y /Adx ~ 1.5 
with Ad x = 0.082i. With increasing disorder, the larger 
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FIG. 4. Normalized s-wave component of the order param- 
eter versus impurity concentration for orthorhombic struc- 
ture. Ago is the magnitude of the order parameter in the 
absence of impurities. 



one, i.e. (Ady), is suppressed faster. When we reach the 
critical impurity concentration both components vanish 
simultaneously. Moreover, <i-wave components of the OP 
decrease linearly with rn and vanish at n% ~ 0.1 as in the 
case of pure d-wave symmetry. Here (■ • ■) denotes aver- 
aging over 20 different impurity configurations. 

As can be seen from Fig. 4, s-wave component also de- 
creases with rii, however while d-wave components ex- 
hibit a linear dependence on impurity concentration s- 
wave component shows a downward curvature similar to 
prediction of the AG theory. 

In conclusion, we investigated the effects of nonmag- 
netic impurities and point defects within a BCS mean- 
field framework by means of BdG equations. For tetrago- 
nal structure, we found out that the observed suppression 
of superconductivity, when impurities are substituted or 
point defects are introduced, can be explained only if the 
OP is d a; 2_j / 2-wave symmetric. In case of s-wave symme- 
try, superconductivity is almost unaffected by disorder. 
When a slight anisotropy is introduced by distorting cop- 
per oxide planes from a square to rectangular lattice we 
observed that a small amount of s-wave contribution is 
induced. For both tetragonal and orthorhombic struc- 
tures we evaluate the critical concentration at which su- 
perconductivity ceases to be very near to experimental 
value ~ 0.1. 
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